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ABSTRACT 

Recent cosmic microwave background anisotropy results from the Wilkinson Microwave Anisotropy 
Probe suggest that the universe was reionized at a redshift around 20 with an optical depth for Thomson- 
scattering of 0.17 ± 0.04. Such an early reionization could arise through the ionizing radiation emitted 
by metal-free population III stars at redshifts of 10 and higher. We discuss infrared background (IRB) 
surface brightness spatial fluctuations from such a generation of early star formation. We show that the 
spatial clustering of these stars at tens of arcminute scales generates a contribution to the angular power 
spectrum of the IRB anisotropies at the same angular scales. This excess can be potentially detected 
when resolved foreground galaxies out to a redshift of a few is removed from the clustering analysis. We 
do not expect faint galaxies at redshifts of ~ 3, with magnitudes less than 20 in the K-band, to be a 
source of strong confusion, since the fractional contribution to the IR background from these galaxies is 
at a level less than a few percent, while the expected contribution from first stars can be more than 50% 
or more. Additionally, assuming a population III stellar spectrum, we suggest that the clustering excess 
related to the first generation of stars can be separated from brightness fluctuations resulting from other 
foreground sources and galaxies using multifrequency observations in the wavelength range of ~ 1 to 5 
jum. In addition to removing the low-redshift galaxy population, the multifrequency data are essential to 
account for certain foreground contaminants such zodiacal light, which, if varying spatially over degree 
scales, can be a significant source of confusion for the proposed study. Using various instruments, we 
study the extent to which spatial fluctuations of the IRB can be studied in the near future. 

Subject headings: cosmology: theory — large scale structure of universe — diffuse radiation — infrared: 
galaxies 



1. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
has now provided strong evidence for an optical depth for 
electron scattering of 0.17 ± 0.04 based on the large scale 
polarization pattern related to rescattering of Cosmic Mi- 
crowave Background (CMB) photons (Kogut et al. 2003). 
If the reionization process is described as instantaneous 
and homogeneous, the measured optical depth implies a 
reionization redshift of ~ 17 ± 5 in a spatially flat uni- 
verse. Such a reionization redshift is higher than previ- 
ously suggested by observational data involving the pres- 
ence of a Gunn-Peterson trough in z ~ 6 quasars (Fan et 
al. 2002). The derived redshift for reionization is at the 
high end of expectations related to reionization scenarios 
based on the ionizing radiation from the first generation 
of star formation associated with population III (hereafter 
Pop III) stars (e.g., Cen 2003; Fukugita & Kawasaki 2003; 
Venkatesan, Tumlinson & Shull 2003; Wyithe & Loeb 
2003). These stars are expected to be very massive with a 
top-heavy mass function (Bromm, Coppi & Larson 1999, 
2002; Abel, Bryan & Norman 2000, 2002), and their de- 
tection is important to understand the astrophysics asso- 
ciated with the reionization process that followed the end 
of dark ages. 

While the direct detection of an individual Pop III star 
or a star cluster is beyond the capability of upcoming tele- 
scopes or instruments, the first star signature can be po- 
tentially detected via indirect methods such as associated 
line emission, e.g., Hell recombination lines of the sur- 



rounding ionized halo (Oh, Haiman & Rees 2001; Tum- 
linson, Giroux & Shull 2001) or through Ly-a emission 
lines (e.g., Tumlinson, Shull & Venkatesan 2003). In Oh, 
Cooray & Kamionkowski (2003), it was suggested that the 
heating of the surrounding interstellar medium by Pop III 
supernovae, and the subsequent transfer of energy to the 
CMB through Compton cooling, lead to a substantial con- 
tribution to small angular scale CMB anisotropies, simi- 
lar to the Sunyaev-Zel'dovich (SZ; Sunyaev & Zel'dovich 
1980) effect. While an individual Pop III supernova re- 
mains unresolved and undetectable, the large bias factors 
of halos containing these supernovae, with respect to the 
linear density field, lead to an excess clustering signature 
in CMB anisotropies at tens of arcminute scales when com- 
pared to the shot-noise contribution related to their finite 
density. The CMB anisotropies related to Pop III super- 
novae can be separated from the dominant SZ fluctuations 
related to galaxy clusters since massive clusters that domi- 
nate the angular power spectrum at arcminute scale can be 
identified and removed in future higher angular resolution 
CMB data (Cooray et al. 2003). 

Here, we suggest that a signature of stars themselves 
can be found in the infrared background (IRB). In par- 
ticular, if Pop III stars are found primarily at redshifts 
between 10 and 30, they are expected to contribute to the 
IRB at wavelengths between 1 and 5 fim. Recent esti- 
mates now suggest that a large fraction of the IRB total 
intensity may be due to these stars (Santos, Bromm & 
Kamionkowski 2002; Salvaterra & Ferrara 2003); A sub- 
stantial IRB arises from the Pop III population not only 
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due to the direct emission associated with these stars, but 
also due to indirect processes that lead to free-free and 
Lyman-alpha emission from the ionized nebulae, or HII 
regions, surrounding these stars. The fractional contri- 
bution to the absolute, and isotropic, background by the 
Pop III population, however, cannot be easily determined 
observationally due to additional contributions from fore- 
ground galaxies, stars, and confusion from zodiacal light. 

Instead of the background intensity, or the monopole, we 
suggest that the presence of first stars can be established 
through spatial fluctuations in the surface brightness of 
the IRB. For this purpose, we make use of the angular 
power spectrum of IRB anisotropics and show that the 
power spectrum is expected to contain an excess cluster- 
ing signature at tens of arcminute scales related to the Pop 
III population. While Pop III stars may dominate tens of 
arcminute-scale fluctuations at short IR wavelengths (be- 
tween 1 and 3 /im), any unresolved galaxies at redshifts 
of 3 and higher, are expected to dominate IRB anisotropy 
at wavelengths greater than 4 /im over all angular scales 
of interest. Thus, we expect the Pop III clustering signa- 
ture can be separated from the angular power spectrum 
due to foreground galaxies and other sources in multifre- 
quency near-IR (NIR) images due to its unique spatial and 
spectral signature. Furthermore, when the angular power 
spectrum due to Pop III stars is estimated, one can use 
that information, for example, to measure the formation 
rate of the first generation of stars in the universe at red- 
shifts between 10 and 30. 

Note that spatial fluctuations in the IRB have already 
been detected in several ground-, suborbital- and space- 
based experiments: Xu ct al. (2002) at 4 /im using a 
rocket-borne experiment; at J, H and K-band with 2MASS 
data by Kashlinsky et al. (2002; also Odenwald et al. 
2003); Infra-Red Telescope in Space (IRTS) data by Mat- 
sumoto 2000; and at large scales with Diffuse Infrared 
Background Experiment (DIRBE) on Cosmic Background 
Explorer (COBE) by Kashlinsky, Mather & Odenwald 
(1996; also, Kaslinsky et al. 1996). Absolute background 
measurements, generally indicate an excess above that ex- 
pected from galaxies alone (Kashlinsky & Odenwald 2000; 
Madau & Pozzetti 2000; Cambresy et al. 2001; Wright & 
Johnson 2001). Additionally, IRTS data also indicate a 
clustering excess at 100 arcmin scales between 1.4 to 2.1 
/im, which remains unexplained by the known properties 
of low redshift galaxies (Matsumoto 2000,2001). Expecta- 
tions for an overall excess in anisotropy include both faint 
galaxies at high redshifts (Jimenez & Kashlinsky 1999) as 
well as a Pop III component related to the first generation 
of stars at redshifts greater than 10 (e.g., Magliocchetti, 
Salvaterra & Ferrara 2003). While simple estimates sug- 
gest that the IRB total intensity can be easily explained 
with Pop III stars alone (Santos, Bromm & Kamionkowski 
2003; Salvaterra & Ferrara 2003), it is useful to also con- 
sider if these stars can be used to explain anisotropy fluc- 
tuations already detected or, if not, the expected level of 
spatial fluctuations from these stars. 

The current observational results related to IRB fluc- 
tuations are, unfortunately, limited to either small angu- 
lar scales (such as ~ 1 to 30 arcseconds measurements 
by 2MASS) or large angular scales (> 0.5 degrees with 
COBE DIRBE), with limited measurements at arcminute 
to tens of arcminute scales. While 2MASS measurements 



by Kashlinsky et al. (2003) are useful to understand the 
presence of a Pop III contribution, at such small angu- 
lar scales, both foreground sources and Pop III stars pro- 
duce power-law like contributions to the angular power 
spectrum and cannot easily be separated (Magliocchetti, 
Salvaterra & Ferrara 2003). Additionally, the Pop III con- 
tribution can also be misidentified easily as a source of 
noise, since, at these small angular scales, one expects a 
shot-noise type power spectrum due to the finite number 
density of these stars. 

In order to understand the clustering nature of the Pop 
III population, here we study the angular power spectrum 
of IRB fluctuations, as a function of wavelength, and intro- 
duce the use of power spectra at each of these wavelengths 
and the cross-spectra between different wavelengths as a 
basis to study the star formation history related to Pop III 
objects. Our calculations extend those of Magliocchetti, 
Salvaterra & Ferrara (2003), by considering the angular 
power spectrum of Pop III stars over a wide range of an- 
gular scales instead of arcsecond scales studied there with 
respect to 2MASS measurements. We argue that a well- 
planned observational program, concentrating on the clus- 
tering at arcminute to tens of arcminute scales or more 
in the wavelength range of 1 to 5 /an, is clearly needed 
to understand the presence of Pop III stars in IRB data. 
Since Pop III stars are expected at redshifts below 30, 
we find that the upcoming and planned wide-field IR and 
near-IR missions from space, such as SIRTF 1 and the pro- 
posed WISE mission, are unlikely to play a major role 
since the Pop III signature is reduced for A > 3 /im. The 
upcoming ASTRO-F 2 mission has imaging capabilities at 
arcsecond resolution in the K, L and M-bands between ~ 
1.7 to 5 microns with its Near Infra-Red Camera (N-IRC) 
and at a higher sensitivity at short wavelengths (Watarai 
et al. 2000; Pearson et al. 2001). While observations 
with ASTRO-F are only restricted to the 10' x 10' field of 
view, these observations are still useful to understand the 
presence of Pop III sources and their contribution to IRB 
spatial fluctuations at angular scales of few arcminutes. 

In order to understand the spatial fluctuations due to 
Pop III sources at tens of arcminute scales and above, and 
also to probe clustering at arcminute scales at wavelengths 
below 2 microns, we also consider a multiwavelength wide- 
field survey using a rocket-borne experiment. This experi- 
ment is planned to image the sky in multiple bands in the 
wavelength range between 1 and 5 /im, with the primary 
goal of exploring the presence of Pop III stars in IRB spa- 
tial fluctuations. While we present a general discussion, 
we will focus on the ability of these experiments to ex- 
tract Pop III information through spatial clustering and 
to remove confusing point sources related to both galactic 
stars and foreground galaxies. 

The paper is organized as follows. In the next section, 
we calculate the power spectrum of IRB anisotropies both 
for Pop III stars and foreground sources. To describe the 
emission from Pop III stars, we make use of the stellar 
spectra calculated by Santos, Bromm & Kamionkowski 
(2003). These calculations also include nebular emission 
associated with both free-free and Lyman-a radiation as- 
sociated from the ionized HII region surrounding individ- 

1 http://sirtf.caltech.cdu/ 
2 http://www.ir.isas.ac.jp/ASTRO-F 
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Fig. 1. — The emission spectrum of a 300 M sun Pop III star at 
a redshift of 15. In addition to the stellar spectrum (solid line), we 
also show the nebular emission related to Lyman-a (dashed lines) 
and free-free (dotted lines) of the ionized HII region surrounding 
the star. This nebular emission depends strongly on the escape 
fraction of ionizing photons from the nebula, / eS c, and we consider 
two extreme cases with / esc = (top lines) and 1 (bottom lines), 
respectively. In the redshift range of 10 to 30, the spectrum peaks 
at wavelengths of order 1 to 3 /im, where the study of Pop III stars 
is strongly favored. 



ual Pop III stars. Note that though we consider Pop III 
sources as the source or ionization, our calculations does 
not necessarily depend on this assumption. The WMAP 
results indicate evidence for a highly efficient source of 
starformation at high redshift and could be in the form 
of Pop II to Pop III stars. While we make use of a Pop 
III spectrum to illustrate the clustering in spatial fluctu- 
ations in the IRB, we also expect similar spatial fluctua- 
tions (though with a different amplitude) if the reioniza- 
tion is related to a different stellar spectrum, say asso- 
ciated with mildly metal-enriched stars, instead of pure 
metal-free stars. In such a scenario, one can simply re- 
place the correct stellar spectrum instead of the Pop III 
one used here and repeat our calculations. When illus- 
trating our calculations, we take cosmological parameters 
from the currently favored ACDM cosmology consistent 
with recent WMAP results (Spergel et al. 2003), including 
a normalization for the matter power spectrum, at scales 
of 8 h~ l Mpc, of cr 8 = 0.9. 

2. ANGULAR POWER SPECTRUM OF IRB FLUCTUATIONS 

Following standard approaches in the literature, we cal- 
culate the angular power spectrum of IRB surface bright- 
ness fluctuations both due to foreground sources at red- 
shifts between and 7 and Pop III stars at redshifts be- 
tween 10 and 30; except for the lower end of the Pop 
III distribution (which we simply take here as the tran- 
sition from Pop III to Pop II stars, and discussed further 
in Section 2.1), our calculations are independent of the end 
points of the redshift ranges considered. 

In general, the clustering contribution to the angu- 
lar power spectrum resulting from foreground galaxies is 



mostly a power-law, which is now well explained with tech- 
niques such as the halo model (e.g., Cooray & Sheth 2002). 
Pop III stars, or any other source that trace linear clus- 
tering at redshifts of order 15 or so, on the other hand, 
have an angular power spectrum that peaks at tens of ar- 
cminute scales. Such a spectrum is expected since high 
redshift sources are expected to be strongly biased with 
respect to the linear dark matter density field and linear 
fluctuation power spectrum contains a peak, or a turnover, 
at the scale corresonding to the matter-radiation equality. 
A detection of such a peak in the power spectrum of any 
tracer source of linear density field would be significant 
since it will help establish our basic uncerstanding of clus- 
tering evolution. Moving to smaller scales corresponding 
to few tens of arcseconds and below, one expects the back- 
ground Pop III stars to show a shot-noise type spectrum 
associated with the finite number density of these sources 
on the sky. In the case of Pop III stars, the transition scale 
from linear to shot-noise type clustering is highly model 
dependent, though with observational data, one can use a 
measured transition as an additional constraint on the as- 
trophysics of the Pop III population. Given that shot-noise 
contributions result from both detector noise and other 
confusions, such as galactic foregrounds, we do not expect 
a precise measurement of the shot-noise power associated 
with Pop III stars to be feasible. For similar reasons, it is 
unlikely that arcsecond-scale IRB clustering information 
is useful for a separation of the Pop III contribution from 
foreground galaxies. 

Our calculational approach to model the angular power 
spectrum of the IRB is similar to the one that was used 
in Knox et al. (2001) to understand the clustering in the 
far-infrared background at wavelengths of few hundred mi- 
crons or more due to dusty star burst galaxies at redshifts 
of ~ 3. The contribution to the IRB intensity, at a given 
wavelength and towards a direction n, can be written as a 
product of the mean IRB cmissivity and its fluctuation 



h(n) = 



dz -^-a(z)j\(z) 
az 



1 



dj x (r(z)n, z) 

h{z) 



(i) 



where ]\{z) is the mean emissivity per comoving unit vol- 
ume at wavelength A as a function of redshift z and r is the 
conformal distance or lookback time, from the observer, 
given by 

r[z) = lHWr (2) 

where the expansion rate for adiabatic cold dark matter 
cosmological models with a cosmological constant is 

n m (i + z) 3 + n K {i + z) 2 + n A ] . (3) 



H 2 = H 2 



Here, Ho can be written as the inverse Hubble distance 
today cHq 1 = 2997.9ft~ 1 Mpc. We follow the conventions 
that in units of the critical density 3Hq/8ttG, the contri- 
bution of each component is denoted fi$, i = c for the 
CDM, b for the baryons, A for the cosmological constant. 
We also define the auxiliary quantities fi m = il c + Of, and 
Qk = 1 — J2i ^it which represent the matter density and 
the contribution of spatial curvature to the expansion rate 
respectively. 

The absolute IRB has now been studied both observa- 
tionally and theoretically (e.g., Hauser & Dwek 2001 for 
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Fig. 2. — Left: A summary of past or upcoming K-band imaging surveys related to clustering studies. The plotted bands represent the 
one-sided sky coverage (with the minimum set at the resolution scale and maximum set at the field of view) in degrees (x-axis) and the survey 
depth between the K-band magnitude at which 10 sources per sky area covered are found down to the 5 a detection limit for resolved sources 
(y-axis). To detect Pop III clustering, surveys that populate the lower left hand corner (large sky-area and large depth) are generally preferred, 
though a combination of experiments, such as a wide-field rocket-borne instrument and the ASTRO-F, may provide a first opportunity. Right: 
Average K-band counts from a variety of sources for both galaxies and stars (based on Saracco et al. 2001). The arrows indicate the 5 a 
magnitude limit to which past or upcoming clustering measurements at K-band can be used to remove foreground sources which are resolved 
in the data, including both instrumental noise and confusion (c.f., Lagache, Dole & Puget 2003 for a description on how to calculate the noise 
based on source counts). In the case of DIRBE, while sources are only resolved down to a higher flux, the clustering analysis accounts for 
most stars based on detailed modeling of the Galactic star distributions (see, Kashlinsky et al. 1996 for details). 



a recent review). Here, we focus on the spatial fluctua- 
tions of the background, SI\(h), and consider the angular 
power spectrum of the IRB, which is simply the Legendre 
transform of the two point correlation function, C{6): 



Ci = 2tt / 6 dOC{9)J (W) 



(4) 



Denoting the Fourier transform of <5/A(n) as <5Ia(1), one 
can define the angular power spectrum of the IRB at wave- 
lengths A and A', in the flat-sky approximation, as 



(rJ7 A (l) 5h (!')) = (2^) 2 M1- 



l')C, AA ' 



(5) 



In order to calculate the spatial fluctuations related to 
the emissivity, we assume Sj\(r(z)ii, z)j j\(z) trace fluctu- 
ations in the source density field, S s = 5p s /p s , such that, 
in Fourier space, 



h{z) 



<5 s (k, z) 



(6) 



The source density field fluctuations are defined by the 
three dimensional power spectrum, which we define as 

(<S.(k, z)<S s (k', z)) = (27r) 3 Mk + k')P ss (fc, z) . (7) 

We will calculate this power spectrum of the source dis- 
tribution, related to both Pop III stars and foreground 
galaxies, based on the halo approach (see, Cooray & Sheth 
2002 for a recent review). In general, the total power spec- 
trum can be written as a combination of 1- and 2-halo 
terms with 
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pih 



ill!' n(m] 



(N s (N s -l)(m)) 



,(k\m) % 



(8) 



where u(k\m) is the density profile (e.g., NFW profile of 
Navarro, Frenk & White 1996) in Fourier space normal- 
ized appropriately with mass, n(m) is the mass function 
(e.g., PS mass function of Press & Schechter 1974), and 
b(m) is the halo bias (e.g., Mo & White 1996; Mo, Jing 
& White 1997). The source distribution within halos is 
encoded by (7V s (m)), the mean source occupation number 
or the average number of individual sources of interest in 
each dark matter halo of mass m, and (N S (N S — l)(m)), the 
second moment of the source distribution. In the case of 
a Poisson-type distribution of sources, (N S (N S — l)(m)) = 
(N s (m)) 2 . In Eq. 8, P hu (k) is the linear power spectrum 
of the density field and we use the transfer function of 
Eisenstein & Hu (1998) to describe the small scale behav- 
ior of this power spectrum. Note that the mean density of 
sources is given by 



n s = / dmn(m)(N s (m)} 



(9) 



Note that this approach to describe clustering of sources 
that lead to the IRB is similar to the one used in Song et 
al. (2003) to describe clustering properties of far-infrared 
background sources. 

Using the Limber approximation (Limber 1954), the an- 
gular power spectrum for a distribution of sources that 
trace a three-dimensional power spectrum P ss (k), when 
projected on the sky, is given by 



dz^^j x (z)j x ,(z)P ss (k = ±.,z) • (101 
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Fig. 3. — The angular power spectrum of IRB anisotropics due to background Pop III stars (top two plots), and background galaxies 
(lower plots). The solid, dashed and dot-dashed lines show clustering at 1 fim, 2 /im, and 4 fim, respectively. In the case of Pop III stars, 
the clustering power spectra are shown for two estimates which are at the high (model A; top-left panel) and low (model B; top-right panel) 
end of the amplitude estimates based on the redshift distribution and biasing factor of dark matter halos, containing Pop III sources, with 
respect to the linear density field. For each of these models, we also plot spectra for the two values of / CS c, with solid lines showing / csc = 
and faint lines showing / csc = 1, to illustrate wavelength dependence on the clustering signal. In addition to the clustering part (2-halo part 
of the power spectrum), we also estimate the shot-noise contribution related to the finite density of Pop III stars. Note that this shot-noise 
contribution related to Pop III stars are highly uncertain due to reasons discussed in the text and we show the upper and lower ranges of 
allowed estimates, each in case of models A and B, respectively, in the top two plots. As shown in the bottom left plot, the foreground galaxies 
have an almost power-law type clustering spectrum; here we only consider galaxies with redshifts greater than 1 under the assumption that 
nearby foreground population can be resolved and removed from the data. Note that we do not show the galaxy shot-noise spectrum here, 
but, refer the reader to Fig. 4 for details. We compare galaxy and Pop III clustering using model A and / eS c = in the bottom-right plot. As 
shown, background Pop III stars are expected to show an excess clustering at tens of arcminute scales due to the strong biasing of dark matter 
halos that contain these stars at redshifts of order 15. In the bottom right plot, for comparison with thick long-dashed lines, we also show the 
total (1-halo + 2-halo) galaxy power spectrum at 2 fim when all galaxies (top curve) and galaxies at redshifts greater than 2 (bottom curve) 
are not removed from the imaging data. The extent to which Pop III clustering can be studied depends on the ability of an experiment to 
resolve galaxies at redshifts out at most to 2 and remove them from the clustering analysis. 



where the comoving angular diameter distance is 

d A = H^VTk 12 smh(H n]l 2 r) . (11) 

Note that as VLk — * 0, <Ia — * r. At I = 50, the Limber 
approximation is valid to within ~ 5% of the exact cal- 
culation that involve radial integrations over the spherical 
Bcssel functions (see, Knox et al. 2001 for the formu- 
lae that involve the exact calculation), and rapidly con- 
verges to the approximated calculation here as one moves 
to higher I. At multipoles below I ~ 50, the Limber ap- 
proximation results in an overestimate of power at the level 
of 10%, and we do not consider this to be important for 
the present discussion given the large uncertainties, which 
lead to several orders of magnitude change in Cj, associ- 
ated with the description of Pop III sources and discussed 
in Section 2.1. 



In addition to the clustering signal, at small angu- 
lar scales, the finite density of sources leads to a shot- 
noise type power spectrum in the IRB spatial fluctuations. 
This shot noise can be estimated through number counts, 
dN/dS, of the contributing sources, as a function of flux 
S, and can be written as 

/■Scut JI\f 

Cf I S M dS > (12) 

where S cu t is the flux cut off value related to the removal 
of resolved sources. 

We will now describe how j\ (z) and ingredients related 
to P ss (k) can be obtained from simple analytical methods 
making use of various approaches in the literature. 
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Fig. 4. — The IRB angular clustering power spectrum at a wavelength of 2 fim. We show our predictions at this wavelength based on 
foreground galaxies (with all galaxies above a redshift of 1) and Pop III stars. The distinct signature of the Pop III population is at tens 
of arcminute scales related to linear clustering at redshifts greater than 10, though its exact amplitude is currently uncertain by at least 
five orders of magnitude based on the optimistic (A) and pessimistic (B) models we have considered. The dashed line shows an estimate 
of the non-linear correction to the Pop III angular clustering power spectrum; non-linearities are important at small angular scales, though, 
at such scales, we expect the shot-noise do be the dominant source of clustering (see, Fig. 3). For reference, we show an average power 
spectrum at 2 (im by Kashlinsky et al. (2003) using deep 2MASS images. The small shaded region surrounding 2MASS data indicate the 
field-to-ficld standard deviation related to the small size of the patch used for these measurements. Following K-band galaxy counts (Fig. 2), 
and the description by Lagache, Dole & Puget (2003), we have calculated the shot- noise contribution related to unresolved galaxies below 
the instrumental sensitivity and show this with a dot-dashed line. An addition contribution beyond this shot-noise may be needed to explain 
2MASS measurements fully, and this could easily involve Pop III sources as was demonstrated by Magliocchetti, Salvaterra &; Ferrara (2003). 
The large hatched region in the upper right corner of the plot show the potential exclusion region for a shot-noise contribution based on 
a representative upper limit on the fluctuations in deep K-band galaxy counts (taken to be K > 24 at 5-<r in a one-arcsecond pixel). Any 
reasonable Pop III clustering contribution, with a corresponding Pop III shot-noise power spectrum that passes through this region, may 
already be constrained. For reference, we also plot the instrumental noise power spectra in K-band of the planned rocket-borne experiment 
and ASTRO-F with instrumental characteristics tabulated in Table 1 and 2; the thick part of these noise curves illustrate the multipoles that 
will be covered due to restrictions from the field-of-view. This rocket experiment will be shot-noise limited by the galactic star and foreground 
galaxy counts and we show an estimate, after removing point sources at the level of 5 a and above, with another dashed-line labeled "SN". 
With 5 <t source removal down to a K-band magnitude of 21.3, the stellar and foreground galaxy shot-noise contribution to the ASTRO-F 
data (shown as a dashed line labeled "ASTRO-F SN" ) is significantly below the expectation for the rocket-borne experiment and will probe 
the presence of Pop III clustering down to the low end of the allowed range in models we have considered. 



2.1. Pop III stars 

To calculate j\{z), we follow the calculation in San- 
tos, Bromm & Kamionkowski (2003) and assume that the 
mass of stars formed per halo is to* = ?7f2f,/f2 m TO when 
m > TO m in and zero otherwise where rj is the star-formation 
efficiency, which we take here as a free parameter. Further- 



more, we assume that Pop III stars trace the star forma- 
tion history at high redshifts and, for calculational pur- 
poses, make use of a model that include molecular hydro- 
gen cooling at a temperature of 400 K. The star formation 
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rate is calculated as 
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where M min is given in Eq. 4 of Santos, Bromm & 
Kamionkowski (2003). The star-formation rate is plotted 
in Fig. 7, and we will return to this later in the context of 
its estimation from Pop III clustering data in the the IRB 
spatial fluctuations. 

Given the star formation rate, we can write 



h{z) 



47T 



(14) 



where (tp op in) = 10 6 years is the mean life time of a 
Pop III star and F\ is the stellar flux as a function of 
wavelength, including nebular emission, of a Pop III star 
(Fig. 1). Again, we make use of calculations by Santos, 
Bromm & Kamionkowski (2003) on the Pop III spectrum 
and consider two extreme cases with / osc = and 1, where 
/esc is the escape fraction of ionizing radiation to general 
intcrgalactic medium from the nebula. The spectra are re- 
produced in Fig. 1 for a star with mass 300 Mq. While the 
stellar spectrum can easily be described by a simple black- 
body, the nebular emission, related to Lyman-a radiation 
and free-free emission, involves a detailed calculation and 
is model dependent on / osc . In the case of / osc = 1, one 
finds more emission at shorter wavelengths and we will see 
later that this leads to an increase in the clustering ampli- 
tude at shortest wavelengths when compared to the case 
with / csc = 0. 

To calculate the three-dimensional power spectrum of 
sources, we make use of the same assumptions as above 
including the fact that the Pop III occupation number is 
simply determined by the halo mass such that (N s (m)) cx 
m when m > m m ; n and zero otherwise. Since the number 
of stars formed is linearly proportional to the halo mass, 
the Pop III clustering is expected to simply trace that of 
dark matter, but with a bias factor determined by the halo 
masses in which stars are found. We also assume a second 
moment which is given by the square of the mean, as in 
the case of a Poisson distribution. At large scales, when 
u{k\m) — ► 1, the simple dependence on the halo mass leads 
to the well known mass averaged bias as 



(brn) = 



jjLt dm m b(m)n(m) 

JZ,. t dm m n ( m ) 



(15) 



where b{m) is the halo bias with respect to the density 
field. 

Since our model has a large number of uncertainties, 
both in terms of the redshift distribution of Pop III 
sources, their spectra, and the biasing factor of halos con- 
taining these stars, we consider a range of models to estab- 
lish both optimistic (model A) and pessimistic (model B) 
estimates of the angular power spectrum. This is similar to 
the approach considered in Oh, Cooray & Kamionkowski 
(2003) to estimate the Pop III supernovae contribution to 
CMB anisotropy fluctuations. Our optimistic estimate in- 
volves: (a) highly biased sources with a cut-off mass for 
Pop III containing dark matter halos that corresponds to a 
temperature of 10 5 K, (b) a low end for the Pop III redshift 



distribution at a redshift of 10, (c) and a star-formation 
efficiency of 100% (r/ = 1) such that all baryons in these 
halos convert to Pop III stars. Our pessimistic estimate 
involves: (a) sources which have low bias factors (with a 
minimum temperature set at 5000 K), (b) the low end of 
the redshift distribution of Pop III sources set at 15, and 
(c) a lower star- formation efficiency with rj = 0.1. In gen- 
eral, the supernovae related to Pop III stars of model A, 
following Oh, Cooray & Kamionkowski (2003), can gen- 
erate the excess CMB anisotropy fluctuations detected by 
small angular-scale CMB experiments. In each of these 
two models, due to the uncertainty related to the Pop III 
spectrum, we also vary / OS c between and 1 to consider 
variations allowed for the amplitude of the angular power 
spectrum. 

For Pop III sources, note that the source power spec- 
trum, in the angular scales of interest, is fully described 
by the linear clustering power spectrum scaled by a bias 
factor. The one-halo term, related to the non-linear part of 
the power spectrum, is not significant due to the fact that 
strong biasing of halos containing Pop III sources leads to 
a dominant two-halo term. This is also consistent with 
the rapid reduction in non-linear clustering one expects at 
rcdshifts of order 10 and higher when compared to today. 
In order to study the extent to which non-linearities may 
be significant, one can introduce the fully non- linear dark 
matter power spectrum based on numerical fits to dark 
matter simulations (e.g., Peacock & Dodds 1996). Since 
sources are not expected to trace dark matter, especially 
in the non-linear regime of dark matter clustering, an al- 
ternative approach is to calculate the non-linear correction 
to source clustering based on the one-halo term related to 
the halo approach. This requires knowledge related to the 
second moment of the Pop III occupation number. Here, 
as an approximation, we make use of the Poisson model 
and assume that the second moment is simply given by 
the square of the first moment. 

In addition to the non-linear clustering, at these small 
angular scales, what determines the power spectrum is the 
shot-noise associated with the finite density of Pop III 
sources. To calculate the shot-noise power spectrum re- 
liably we need a detailed model describing number counts 
of Pop III sources, including the exact mass function of 
Pop III stars within each of the dark matter halos. For 
the clustering calculation, this information is not required 
since, in the end, the clustering amplitude depends only 
on is the total baryonic mass within each halo converted 
to stars, and not how this mass gets divided in to stars. 
Additionally, the shot-noise contribution is more uncer- 
tain than the clustered component since the shot-noise de- 
pends on the rarest Pop III sources which also happen to 
be the brightest. Similar to our approach for the cluster- 
ing signal, we also make both optimistic and pessimistic 
estimates of the shot-noise contribution, but for simplic- 
ity, we only show the shot-noise related to the / CS c = 
spectrum (see, Fig. 3); in cases where the clustering sig- 
nal is high due to higher biasing of rarer halos, we also 
find a higher shot-noise contribution since such rare ha- 
los have a low projected surface density. Note that our 
highest shot-noise estimates, corresponding to the most 
optimistic estimate on the Pop III clustering contribution, 
can be excluded by fluctuation analyses of already exist- 
ing deep K-band imaging data. Thus, the recent 2MASS 
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Fig. 5. — The correlation coefficient of IRB power spectra. Here, 
we show the cross-correlation between 1 (im map and the maps at 
2, 4 and 8 fim, under the assumption that Pop III sources follow 
our most optimistic description (model A) and that resolved fore- 
ground galaxies are removed out to a redshift of unity. The top 
lines show the correlation if galaxies were the only contributor to 
IRB anisotropics. The Pop III clustering only lead to a correlation 
at tens of arcminute scales, which when the two contributions are 
added, result in an overall reduction of the total correlation coeffi- 
cient at arcminute scales. 

clustering measurement at small angular scales, discussed 
with respect to figure 4, cannot be explained with Pop III 
shot-noise alone. 

2.2. Foreground galaxies 

To describe the spectrum of an individual galaxy, we 
make use of the spectrum of M82 (e.g., Silva et al. 1998 
where a model fit is presented to explain the observed spec- 
tral energy distribution over a wide range of wavelengths 
from far- infrared to UV light). While the use of such a 
spectrum to describe all galaxies is likely to be an overes- 
timate, we consider it here since this provides a conserva- 
tive upper limit on the extent to which foreground galaxies 
can confuse clustering studies of the first star population 
using IRB data. To describe the redshift distribution of 
foreground galaxies, we make use of the normalized star- 
formation rate between redshifts of and 7, but intro- 
duce an exponential cut off at redshifts above 4 such that 
the redshift distribution turns over and converges essen- 
tially to zero at a value of 7. Note that the observational 
data related to galaxy counts, at redshifts greater than 3, 
are highly uncertain, though preliminary studies indicate 
that the fractional contribution to the infrared background 
from galaxies at redshifts greater than 3 (and K-band mag- 
nitudes fainter than ~ 23), may not be significant (Totani 
et al. 2001a). Thus, we expect our modeling to provide 
only an estimate of the galaxy contribution to the infrared 
background. While a lack of massive halos at redshifts 
greater than 3 reduces the presence of a large number of 
IR-bright galaxies considerably, the presence of unusually 
red galaxies in certain K-band imaging data (e.g., Totani 
et al. 2001b; Daddi et al. 2003) suggests that the situation 
may be more complex than suggested by numerical calcu- 



lations. In addition to the suggested distribution, we also 
considered alternative redshift distributions for foreground 
galaxies, but found consistent results as long as the distri- 
bution is not sharply peaked towards the low redshift end, 
such as below z of 1; such models are clearly inconsistent 
with K-band galaxy counts among other data and we do 
not consider such models further. 

The clustering of foreground galaxies is also calculated 
following the halo-based approach with a halo occupation 
number based on fits to semi-analytic numerical simula- 
tions data by Sheth & Diaferio (2001): 

(N) _fN lQ ll M & h- l <m<M B , . 

^ s > ~ \ N (m/M B ) a m > M B [W) 

where we set Mb = 4 x 10 12 Mq/H and take a and Nq as 
free parameters. In our fiducial description, these two pa- 
rameters take numerical values of 0.8 and 0.7, respectively. 
In addition to the mean, we also account for departures 
from the mean in the 1-halo term with a detailed model for 
the second moment, instead of the simple Poisson descrip- 
tion used in above. Following Scoccimarro et al. (2001), 
we make use of the binomial distribution, matched to nu- 
merical data, to obtain a convenient approximation. The 
second moment is then 

(N s (N s -l)y/ 2 =(3(m)(N s ), (17) 

where /3(m) = log y/m/W 11 ^ 1 M Q for m < 1O 13 /i -1 M 
and (3(m) — 1 thereafter. Note that in equation 8, the sim- 
plest approach is to set p = 2 when calculating P^(k). 
In halos which contain only a single galaxy, however, we 
assume that this galaxy sits at the center and take p = 2 
when (N S (N S — 1)) is greater than unity and p = 1 other- 
wise. These description for galaxy clustering lead to a 
power-law like three dimensional spectrum over a wide 
range of physical scales and also explains why the clus- 
tering of galaxies shows no excess, as in dark matter clus- 
tering, between the transition from linear to non-linear 
clustering (see, Cooray & Sheth 2002 for details). 

While the use of a low-redshift determined occupation 
number, motivated by semi-analytical modeling of "blue" 
galaxies, to describe clustering properties of IR galaxies 
may seem problematic, we suggest that this description 
is adequate for purposes of the present discussion. The 
main problem related to our description comes from the 
fact that, based on observations, one naively expects the 
minimum mass at which an IR galaxy at redshifts of ~ 3 
to form in halos with a mass at least two orders of mag- 
nitude above the minimum mass we have indicated in Eq. 
16. For example, the existing clustering data of IR galaxies 
at high redshifts (e.g., Daddi et al. 2003) indicate strong 
clustering of K-band selected galaxies at a level similar to 
or greater than the clustering strength of the Lyman break 
galaxy population at redshifts ~ 3 and suggests that the 
minimum mass must be set at a level of ~ 10 13 Mq. In the 
case of Lyman-break galaxies, the halo occupation number 
is investigated in detail by Bullock, Wechsler & Somerville 
(2002). Contrary to naive expectations, however, these au- 
thors find an occupation number with a minimum mass at 
the level we have indicated in Eq. 16, though there is a 
large range of allowed values (see, e.g., Fig. 3 of Bullock, 
Wechsler & Somerville 2002). The difference can be un- 
derstood based on the fact that the minimum mass, as 



9 



indicated by the halo occupation number, is not necessar- 
ily the one related to the mass above which halos are found 
to contain one observed object, but rather the minimum 
mass in which one galaxy may form. While analytically, 
halos with mass below the minimum expected and above 
the minimum allowed by the halo occupation number form 
a fraction of a galaxy in each halo, in reality, this is meant 
to indicate the fact that only a fraction of such halos form 
galaxies occasionally. On the other hand, all halos with 
mass above the expected minimum mass, which in this case 
is two orders of magnitude higher than the inferred min- 
imum, are expected to host one or more galaxies. While 
the overall dilution of the clustering strength by a sample 
of small mass halos, where galaxies are occasionally found, 
is not significant, these halos are important to describe the 
clustering properties of the high redshift population (see, 
Bullock, Wechsler & Somerville 2002 for a further discus- 
sion). 

We note that our predictions related to IR galaxy clus- 
tering at redshifts around 2 are consistent with Daddi 
et al. (2003) measurements and find this to be adquate 
enough for the present paper. This is mainly due to the 
fact that we are mainly concerned with clustering prop- 
erties of brightness galaxies at redshifts between 1 and 2, 
or in the case of K-band, at magnitudes below 20 or so. 
As we will soon discuss, this is within the range of mag- 
nitudes for which upcoming experiments will be able to 
remove resolved sources from a clustering analysis. We 
are less concerned about the clustering properties of faint 
high redshfit galaxies, at z ~ 3, for one significant reason. 
The suggested measurement involves a clustering analy- 
sis of brightness fluctuations and not source counts, as in 
traditional studies of galaxy clustering at low redshifts. 
Thus, the clustering power spectrum from each redshift 
range, or bin, is weighted by the contribution to the abso- 
lute IR background produced by galaxies in that range or 
bin. Assuming fluctuations in the brightness trace source 
counts, we can roughly write this as Ci ~ I 2 wi where wi 
is the angular power spectrum of source clustering. The 
fractional contribution to IR background from high red- 
shifts is substantially less than that produced by nearby 
bright galaxies out to a z of 1 (e.g., Cambresy et al. 2001). 
In the K-band, the fractional contribution to the infrared 
background from galaxies with magnitudes fainter than ~ 
20, which are expected to be at redshifts greater than 2, 
is at most at the level of 10%, or if not less. On the other 
hand, the fractional contribution from Pop III sources can 
be at the level of 50%, or more, especially under extreme 
models discussed in the literature (e.g., Salvaterra & Fer- 
rara 2003). In relative terms, using the fractional contri- 
bution to the background, one can safely conclude that 
Pop III clustering will dominate brightness fluctuations, 
when compared to that of faint galaxies at redshifts of 
around 3, by at least a factor of several tens or more since 

c PopIII /c K>20 K [7 PopIn // K 

>2o] • The extent to which 
this ratio can be known exactly depends on the relative 
strength of source clustering at redshifts of around 3 when 
compared to that at a redshift around 15. Though some 
subsamples of faint galaxies, such as the red galaxy (J- 
K > 1.7) population of Daddi et al. (2003) at redshift 
of 3, are strongly clustered, we do not necessarily expect 
this clustering strength to be significantly higher than Pop 



III sources such that C t op < C l , . However, we note 
that this conclusion could substantially change if the z ~ 3 
population is to be a significant contributor to the IR 
background and its fractional contribution is substantially 
higher than 10% or so level implied in the Cambresy et al. 
(2001) paper. 

The shot-noise power spectrum related to galaxies, as 
well as galactic stars, is more certain due to the avail- 
ability of the number counts in certain wavelength ranges 
between 1 and 5 microns. In Fig. 2 (right plot), we plot 
the average number counts in the K-band compiled from a 
variety of sources in the literature following Saracco et al. 
(2001) and Cambresy et al. (2001; and references therein). 
We will use these number counts to establish the confu- 
sion noise in various experiments, summarized in Fig. 2 left 
plot, at near IR wavelengths to understand the extent to 
which these foreground sources dominate clustering stud- 
ies through their shot-noise power spectrum. In both plots 
of Fig. 2, we also illustrate the magnitude limit to which 
resolved sources can be removed at the 5 a instrument 
noise level in various past and upcoming wide-field imag- 
ing data that are adequate enough for a clustering analysis 
of any first-star contribution to the IRB. 

2.3. Angular Power Spectra 

In Fig. 3 (top plot), we show the angular power spectra 
of galaxies and Pop III sources as a function of wavelength 
at 1, 2 and 4 /im. These wavelengths are loosely chosen 
to represent the clustering behavior of sources in the IR 
regime. In the case of galaxies, the angular power spec- 
trum is nearly a power-law over a wide range of angular 
scales. The excess clustering at degree scales is related to 
the linear power spectrum which projects at low redshifts 
(z ~ 1). When interpreting this galaxy power spectrum, 
one should consider this as a representative case where 
foreground galaxies, out to a redshift of ~ 1, is removed 
from the data such that nearby low redshift galaxies are 
absent; In our calculations, we allow for such a removal 
by simply cutting off the redshift distribution of galaxies 
at redshifts below 1. Note that our approach to describe 
clustering of galaxies at IR wavelengths is different from 
previous approaches. Instead of extending the measured 
power spectrum at low redshifts, such as from the APM 
survey, to higher redshifts, or using numerically calibrated 
fitting functions for the non-linear power spectrum, such 
as from Peacock & Dodds (1996), which are valid only for 
dark matter, we have made use of the halo approach to 
calculate the clustering of galaxies directly. 

Our approach to describe clustering is also independent 
of techniques that make use of galaxy number counts at 
IR wavelengths. While there is now detailed information 
related to galaxy counts, the exact redshift distribution of 
IR contributing galaxies still remains somewhat uncertain. 
As a simple approach, here we have made use of a distri- 
bution based on the normalized star-formation rate. As 
long as this distribution is not dominated by galaxies at 
the high end of this redshift distribution, galaxy clustering 
follows a power-law. Though we have not considered here, 
our halo-based approach can be extended, under certain 
assumptions, to model number counts of IR galaxies simi- 
lar to approaches considered in studies such as by Jimenez 
& Kashlinsky (1999). 

The Pop III stars, since they are present at redshifts 



10 



greater than 10 and are highly biased with respect to the 
density field, trace the linear regime of clustering at tens 
of arcminute scales when projected on the sky today. As 
shown in Fig. 3 top plots, the Pop III stars, thus, show the 
expected signature of an excess clustering, above the shot- 
noise power spectrum, at angular scales where the linear 
power spectrum has a peak associated with the matter- 
radiation equality. The difference between the linear clus- 
tering of galaxies at low redshifts and at angular scales of 
few degrees or more, and Pop III sources at high redshifts, 
and at tens of arcminute scales, can be understood based 
on the redshift evolution of the linear density field power 
spectrum alone. As discussed in Cooray et al. (2001), this 
evolution can in fact be used to constrain certain cosmolog- 
ical parameters though we do not pursue such possibilities 
here. 

In Fig. 3, we also show Pop III source clustering angu- 
lar power spectra for flux spectra involving / csc = and 1. 
We consider these two values as a potential bound on the 
uncertainties related to the Pop III spectrum. While the 
clustering excess at a short wavelength like 1 /xm is small 
in the case of / CS c = 0, with / csc = 1, we find a signifi- 
cant contribution at the level comparable to 2 fim. This 
can be understood based on the Pop III spectra we show 
in Fig. 1 and results from the fact that when / csc = 1, 
the short wavelength cut off is decreased to a wavelength 
value lower than for / csc = 0. Additionally, due to the de- 
crease in free- free emission with / csc = 1, one finds a lower 
clustering amplitude at high end of the IR wavelengths 
when compared to the case with / esc = 0. While we have 
considered these two extreme cases one expects / csc to be 
order few tens of percent at most, such that more realis- 
tic scenario will be more closely related to / csc = than 

/csc 1 • 

In addition to the spectrum, other uncertainties in 
physics of the Pop III population lead to a highly indeter- 
minate overall normalization for the angular power spec- 
trum. These uncertainties, which are more significant than 
/csc alone, are related to the Pop III redshift distribution 
and biasing with respect to the linear density field. We 
bound these uncertainties and consider two estimates at 
the high and low end; it is likely that the true cluster- 
ing is somewhere between the two, though, due to large 
uncertainties, we are unable to determine this to a better 
accuracy than the range implied by the two models. While 
the amplitude of the Pop III angular power spectrum may 
be uncertain, however, the angular scale at which the Pop 
III clustering signature is expected is more reliable as it is 
simply a reflection of the linear power spectrum projected 
at redshifts of order 20; varying the projection from red- 
shifts 10 to 50 or so does not lead to a large change in this 
angular scale since the angular diameter distance does not 
change significantly in this range of redshifts. 

In addition to a clustered signal, the finite density lead 
to a Poisson noise. These shot-noise power spectra for Pop 
III sources, as a function of the wavelength, is also shown 
in Fig. 3 (top plots). As we discussed earlier, the shot-noise 
is more uncertain than the clustering power spectrum since 
the shot-noise depends strongly on the rare and bright 
events while the clustering amplitude, as well as the cu- 
mulative background, does not. While the shot-noise can 
easily be confused with the shot-noise power due to faint 
unresolved galaxies, among other foregrounds, the cluster- 



ing excess, at tens of arcminute scales, however, provides 
a potentially interesting signature of Pop III which can be 
observed. 

As shown in Fig. 3 (bottom right plot), the excess Pop 
III clustering signature lies above the angular power spec- 
trum of galaxies when wavelengths are in the range of ~ 1 
to 2 /Ltm and when nearby galaxies are removed from the 
data. Allowing for the presence of nearby galaxies leads 
to a slightly higher clustering signal for galaxies, which we 
have shown with a long-dashed line in the right-bottom 
plot of Fig. 3 for a wavelength of 2 (ira. The removal of 
bright nearby galaxies leads to a decrease in the non-linear 
clustering, but not necessarily the large angular scale clus- 
tering related to the 2-halo part of the power spectrum 
since the large angular scale clustering is determined by 
the linear power spectrum scaled by the bias factor. For 
galaxy models considered here, the average bias factor is 
not strongly sensitive to the removal of nearby sources. In 
general, at wavelengths above a few microns, the galaxy 
clustering fully dominates IRB anisotropies, though, one 
can improve the detection of Pop III stars even at these 
wavelengths if most resolved galaxies are removed. For 
comparison, in Fig. 3, the lower thick long-dashed line 
show the expected clustering in the infrared-background at 
2 /im from galaxies at redshifts greater than 2 suggesting 
that a removal down to this redshift range (or magnitudes 
at the level of 25 in K-band) is essential to extract the 
clustering signal associated with first-star population. In 
general, to perform such a removal and to detect Pop III 
clustering signal at tens of arcminute scales, an experiment 
with subarcminute-scale angular resolution, but covering 
several tens of square degrees, is needed. 

In Fig. 4, we highlight the main aspects of Pop III clus- 
tering. For reference, we also show the 2MASS clustering 
measurement in the K-band at small angular scales (Kash- 
linsky et al. 2003) as a reference. Note that the Pop III 
signature is present at tens of arcminute scales with an 
excess above the level expected from galaxies (assuming 
that foreground galaxies, out to a redshifts <~ 1 to 2, can 
be resolved and removed from the analysis). In addition to 
2MASS, DIRBE (Kashlinsky, Mather & Odenwald 1996) 
and IRTS (Matsumoto 2000) measured Pop III clustering 
at large angular scales and detected excess fluctuations 
(discussed with respect to Fig. 6). In addition to the lin- 
ear clustering of Pop III sources, we also show the correc- 
tion to the power spectrum associated with non-linearities. 
While non-linearities is expected to increase the power by 
an order of magnitude at very small scales (I ~ 10 5 to 10 6 ), 
a comparison to Fig. 3 reveal that at such small angular 
scales, the clustering spectrum related to Pop III sources is 
more likely to be dominated by their shot-noise. At these 
same angular scales, the shot-noise related to unresolved 
foreground galaxies also become important such that the 
individual contributions related to galaxies and first stars 
are not easily separable. 

Based on our description of the Pop III population and 
its clustering aspects, we find the clustering amplitude re- 
sulting from emission by first stars is not fully responsi- 
ble for the clustering signal detected in the 2MASS data 
at arcsecond angular scales. The 2MASS angular power 
spectrum is more likely due to a combination of the non- 
linear clustering of unresolved galaxies and the shot-noise 
associated with such sources. The number counts of galax- 
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ies in the K-band now extend down to a magnitude of 24, 
and if any shot-noise existed at the 2MASS level, its fluc- 
tuation amplitude from pixel to pixel, due to the Poisson 
behavior, would correspond to a magnitude ~ 22 to 23 in 
one-arcsecond pixels. The deep counts will easily be sen- 
sitive to such variations and based on such an argument, 
we suggest that the Pop III source shot-noise or clustering 
contribution must be below the large shaded region on the 
top right corner of Fig. 4. This region essentially corre- 
sponds to 5 o fluctuations down to a K-band magnitude 
of 24 in a one-arcsecond pixel. On the other hand, using 
measured galaxy counts alone, a significant shot-noise is 
expected from unresolved sources that are below the in- 
strumental detection limit but above the confusion limit 
(e.g., Lagache, Dole & Puget 2003 for a description); we 
show this contribution, together with the shot-noise from 
sources below the confusion limit, as a dot-dashed line in 
Fig. 4. It is likely that a combination of this shot-noise 
and the non-linear clustering of the unresolved sources 
contribute to the 2MASS measurements, though, a sub- 
stantial contribution from non-linear clustering of Pop III 
sources may be expected in certain descriptions of the Pop 
III population (Magliocchetti, Salvaterra & Ferrara 2003). 

Since various possibilities may exist to explain 2MASS 
measurements (either in terms of galaxies, Pop III sources 
or in comibation), we suggest that further investigation 
of small scale clustering behavior of the IRB may not be 
adequate to separate these contributions. While one can 
extend a survey such as 2MASS to large angular scales 
by tiling subsequent fields, this can cause problems when 
combing adjacent fields to measure large scale clustering 
beyond that of a single field of view (see discussions in 
Odenwald et al. 2003). The alternative approach is to 
image a wider area (of tens of square degrees) directly at 
slightly lower resolution and perform a clustering analysis, 
either after removing resolved stars and galaxies directly 
in such an image or after accounting for such sources based 
on higher resolution data of smaller fields, to see if there 
is any evidence for an excess as indicated by Fig. 4 related 
to Pop III sources. The best combination to search for 
the Pop III signal appears to be a wide-field camera with 
fidelity on degree scales combined with a deep survey, such 
as ASTRO-F, to minimize foreground confusions. 

The detection of excess Pop III source clustering will be 
limited to the extent that foreground stars and galaxies 
can be removed. In surveys where the source removal is 
not significant, one expects a higher shot-noise contribu- 
tion from galactic stars and foreground galaxies, and this 
shot-noise level could potentially be above the level of ex- 
pected clustering of galaxies. To model the expected shot- 
noise, we follow calculations by Lagache, Dole & Puget 
(2003), and we show the case for a wide-field rocket-borne 
experiment (see Table 2 for details) in Fig. 4. The in- 
strumental noise power spectrum C^ oiso related to each 
experiment is shown with a dashed-line and is given by 

Cr se = ^Uy^e l2{Ae)2 , (18) 

where <7 p i x is the noise-per-pixel (see, Table 1), N p - lx is the 
number of pixels and A9 is the pixel scale or the scale of 
additional smoothing, if the latter is employed. As writ- 
ten, the noise contribution to the measured power spectra 



at each wavelength is described through standard CMB 
analysis approaches (Knox 1995) given the filtered Gaus- 
sian beam and the noise of the whole focal plane array. 
With the noise power spectrum, one can write the error, 
or standard deviation, for a measurement of the angular 
power spectrum at each multiple as 

sc <=\0^[ c < +c ' n °'*]- <i9 » 

where the term related to Ci accounts for the cosmic vari- 
ance associated with the finite sky-coverage. When the 
power spectrum measurements are binned in the multi- 
pole space, the error associated with the binned power 
spectrum measurement is reduced further by a factor of 
VAl where Al is the width of the bin. While we have 
not shown these errors in Fig. 4, due to uncertain C; of 
Pop III stars, the extent to which upcoming instruments 
probe the power spectrum, related to / s k y is shown with a 
combination of thick and thin lines when plotting Cf omc . 

For instrumental parameter values in Table 1 related 
to the rocket-based imaging data, one can remove point 
sources, at the 3 er confidence level above the expected 
instrumental noise and confusion noise from sources below 
the resolution limit, down to a flux level of ~ 2 x 10 -7 nW 
m~ 2 . This corresponds to a K-band magnitude of 16.8. 
Using K-band counts for stars and galaxies, we expect a 
shot-noise at the level of C\ <~ 5 x 10~ 8 nW 2 m~ 4 sr _1 due 
to unresolved and unsubtracted sources. We show this 
contribution with a dashed line labeled "SN" in Fig. 4. 
In the case of a wide-field rocket-borne experiment, one 
expects a potential detection of the Pop III population, at 
multipoles of order 10 3 , if these stars follow our models at 
the optimistic end. 

At the low end of our models, the Pop III clustering 
amplitude is lower and the detection will strongly depend 
on the extent to which galaxies, with K-band magnitudes 
at the level of 20 and fainter, can be resolved and removed 
from the data. While reaching such a low magnitude level 
for resolved sources is beyond the capability of the rocket- 
based experiment, due to restrictions in the integration 
time among others, such a study is clearly possible with 
ASTRO-F. At the 5 a confidence level above instrumental 
noise, ASTRO-F allows identification and removal of fore- 
ground sources down to a magnitude limit of 21.3 in the 
K-band or 2.2 x 10~ 9 nW m~ 2 in 1.4 arcsecond-pixels. 
Using the same counts as before, we now find a shot-noise 
contribution at the level of Ci ~ 2.0 x 10~ 9 nW 2 m" 4 
sr _1 ; while the extent to which Pop III star clustering 
can be detected with ASTRO-F depends on the cluster- 
ing nature of unresolved galaxies below its detection limit, 
we find that ASTRO-F is useful for the study of Pop III 
stars since it allows a large range in the clustering ampli- 
tude to be probed at arcminute scales due to the limited 
ficld-of-view, 10 arcmins by 10 arcmins. Such an angular 
power spectrum measurement is useful to understand the 
extent to which clustering of Pop III sources contributes 
to spatial fluctuations of the IRB given that we predict the 
presence of an excess at angular scales close to 10 arcmins 
and above. In addition to the limited field-of-view, an- 
other drawback with ASTRO-F is the limited wavelength 
coverage restricted to K-band and longer. To properly es- 
tablish the clustering of Pop III sources, one will require 



12 



a study of the frequency spectrum related to the clustered 
component and lack of observations below K-band may 
complicate such an analysis; Furthermore, if coverage at 
shorter wavelengths existed, one can use such data for an 
extraction of the Pop III contribution, based on its ex- 
pected spectrum, such that the confusion between Pop III 
sources and foreground galaxies are further reduced. 

While observations in a single band, e.g., K-band, is 
useful to understand Pop III clustering, this can be better 
achieved with multiwavelength data as different channels 
can be combined to study clustering in each band and be- 
tween bands. In Fig. 5, as an illustration, we show the cor- 
relation coefficient related cross wavelength power spectra. 
This correlation coefficient is defined as 

C{\ i ,\ j ) = ^£==, (20) 

under the assumption that the Pop III sources follow our 
most optimistic assumptions (model A) and that the fore- 
ground galaxies are removed out to a redshift of 1. As 
shown, the total correlation coefficient ranges near unity at 
large angular scales, decreases at tens of arcminute scales, 
and increases back to almost unity at small angular scales. 
If the Pop III contribution were to be non-existent, the cor- 
relation coefficient would be more uniform while the de- 
crease in the correlation at arcminute scales is associated 
with the Pop III population. While an individual angular 
power spectra measured at each wavelength was expected 
to show excess clustering due to Pop III stars, the corre- 
lation coefficient from cross-power spectra is expected to 
show a decrease at the same angular scales. 

In addition to separating foreground galaxies and back- 
ground first-star contribution to spatial clustering in the 
infrared background, multiwavelength observations are de- 
sirable to understand additional contaminants, among 
which zodiacal light is important. While spatial clustering 
properties of zodiacal light, or dust particles responsible 
for this emission especially at degree angular-scales, is un- 
known, one expects a distinct spectral dependence to its 
emission given the size-distribution of dust particles. With 
observations at several or more frequencies in the IR band, 
one can make use of this expected spectral shape of the 
zodiacal emission to partly remove its contribution and to 
preform a clustering study related to the zodiacal emission 
separately. This is similar to what has been done in esti- 
mating the total IR background in imaging data such as 
with DIRBE (e.g., Kashlinsky et al. 1996). With imag- 
ing data such as from a rocket, another approach is to 
make use of the temporal variations related to the zodi- 
acal emission. Since the zodiacal light is associated with 
micron-size particles near the Earth's orbit, observations 
spread over the Earth's orbit that span few months or 
more are expected to view a different dust content to- 
wards the same background sky-area at different obser- 
vational epochs. This could give additional information 
to remove the zodiacal light, though, we emphasize that 
the best option to remove this contaminat is to perform 
multiwavelength imaging data over the IR band. 

As a general consistency check on our models, instead 
of calculating the angular power spectrum, we also calcu- 
lated the mean or the absolute background to be compared 



with observed results. Note that our Pop III redshift dis- 
tribution follows that of Santos, Bromm & Kamionkowski 
(2003) and since we use the same stellar spectrum as they 
calculated, the Pop III star model produces the absolute 
background consistent with Fig. 7 of their paper, to the 
extent that star-formation efficiency is taken here as a free 
parameter, and by varying the minimum redshift at which 
Pop III sources are found to 7 from 10; in our fiducial de- 
scription of the Pop III population, we use a some what 
higher minimum redshift for Pop III sources and the IR 
background related to this model results in an overall de- 
crease in the Pop III contribution at wavelengths below 
~ 1.5 /im. Our use of a higher minimum redshift for Pop 
III sources is motivated by arguments related to the ex- 
pected rapid transition from Pop III to Pop II sources, 
though, we note that based on calculations for Pop II star 
spectra by Tumlinson & Shull (2000), Pop II stars will con- 
tinue to make a substantial contribution to the IR back- 
ground and background fluctuations at lower wavelengths. 
In the case of foreground galaxies, our model produces less 
than ~ 40% of the background at 2 /im; this is consistent, 
though potentially higher than, estimates based on galaxy 
number counts (Cambresy et al. 2001). Note that a sig- 
nificant fraction of this contribution comes from sources at 
rcdshifts below <~ 2 and, as expected from studies involv- 
ing source counts, we expect galaxies out to redshift <~ 2 
to be the dominant confusion for a study such as this and 
the extent to which such galaxies can be removed from the 
data will determine the ability to extract Pop III clustering 
excess in the surface brightness fluctuations. 

In order to consider the extent to which our cluster- 
ing predictions are consistent with prior results, we com- 
pare our angular power spectra based on predictions with 
COBE DIRBE measurements, as a function of wavelength, 
by Kashlinsky, Mather & Odenwald (1996). These au- 
thors presented a measurement of C(6 = 0), the correla- 
tion function at zero lag, or rms fluctuations at the beam 
scale, at several wavelengths and using a top-hat filtered 
maps. This quantity can be calculated from our power 
spectra using the fact that 

C(0) ee (5ll) a = J JlL C ^'w\1<j) , (21) 

where W(x) = 2J\(x)/x is the top-hat window func- 
tion and a — 0.46° is the angular scale of filtering. We 
summarize our results in Fig. 6, where we plot y/C(0) 
for the three cases involving galaxies and Pop III sources 
(with ,/c SC = and 1), and a comparison to COBE DIRBE 
measurements at 1.25, 2.2, 3.5 and 4.9 /im. Note that, in 
general, our model predictions, even at the optimistic end 
for Pop III sources, suggest at least an order of magnitude 
lower values for \/C(fi) than measured. The prediction 
based on galaxies alone is similar and is consistent with 
previous estimates (e.g., Kashlinsky ct al. 1996). While 
we have some freedom to fit rest of the data by varying our 
model parameters, we have not attempted to perform such 
an analysis here. The Pop III sources, if at the low end of 
our prediction, suggest roughly three orders of magnitude 
smaller contribution to clustering seen in DIRBE data. 

For comparison, in Fig. 6, we also show results related to 
the fluctuation analysis in IRTS data as a function of the 
wavelength in the IR regime. These fluctuation measure- 
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Fig. 6. — The correlation function at zero-lag and a comparison to 
DIRBE results from Kashlinsky, Mather & Odenwald (1996), IRTS 
results by Matsumoto ct al. (2000, 2001), and NITE fluctuation 
measurement by Xu et al. (2002); note the difference in angular 
scales at which these measurements are made: DIRBE measure- 
ments are at a filtering scale of 0.46°, IRTS measurements are at 
the beam scale of 8 arcmins by 20 arcmins and the NITE measure- 
ment is at an angular scale of 17 arcseconds. Except at 2 (im, note 
that our model predictions based on galaxies are roughly an order 
of magnitude lower than the measurements. The expected contribu- 
tion from Pop III stars span at least three orders of magnitude and 
arc highly uncertain; here, we show the two ranges implies by our 
data with model A (top line), with / esc = 0, and model B (bottom 
line), with / esc = 1. 



Table 1 



Instrumental Characteristics; ASTRO-F 



Band 


A 


A0 


5 F 


f^pix 




(Mm) 


(arcsec) 


/ijy (5-cr) 


(nW m~ 2 sr- 1 ) 


K 


1.8-2.7 


1.46 


1.3 


14.5 


L 


2.7-3.7 


1.46 


1.2 


15.6 


M 


3.7-5.1 


1.46 


1.9 


17.6 



NOTES. — The sensitivity, S F, and noise-per-pixel, <r p i x , 
is based on 500 seconds of integration per pointing. The 
Near Infra-Red Camera (NIRC) of the ASTRO-F mission 
contains arrays of 512 x 412 pixels with an effective area 
of 10 2 arcmin 2 per pointing. See, Watari et al. (2000) for 
more details. 



Table 2 



Instrumental Characteristics; Rocket-Borne 
Experiment 



Band 


A 


AA/A 


S F 


Cpix 




(H 




AtJy (5-cr) 


(nW m~ 2 sr 


J 


1.25 


0.24 


530 


7.5 


H 


1.65 


0.17 


325 


6.5 


K 


2.2 


0.16 


205 


4.4 


L 


3.5 


0.26 


134 


1.8 



ments are made at the IRTS beam scale of 8 arcmins by 20 
arcmins, which probes slightly smaller angular scales than 
the DIRBE measurements shown in Fig. 6. There is con- 
siderable agreement, however, between DIRBE measure- 
ments and IRTS. It is unclear if these excess fluctuations 
are due to the shot-noise associated with unresolved galac- 
tic stars, due to the large beams of these experiments, or 
an additional extra-galactic component. Note that these 
observations probed large angular scales, while a substan- 
tial contribution from the Pop III population is expected 
at tens of arcminute scales and, thus, these data do not 
necessarily aid in understanding the extent to which first 
stars are an important source of spatial fluctuations in IR 
wavelengths. In future cxpriments attempting to study 
Pop III clustering, sufficient angular resolution to measure 
and remove bright stars and galaxies is essential. 

3. RECONSTRUCTING THE STAR FORMATION HISTORY 

Since the Pop III redshift distribution is expected to 
closely follow the star formation history at redshifts be- 
tween 10 and 30 or so, the anisotropy power spectrum of 
the IRB can be used to extract information on the star for- 
mation history, under the assumption of certain aspects 
related to the source clustering and the spectrum. For 
this purpose, we make use of multi-frequency measure- 
ments including the cross power spectra between frequen- 
cies. To understand the extent to which star formation 
history can be studied with IRB fluctuations, we make 



NOTES. — The sensitivity, S F, and noise-per-pixel, <r p i x , 
assumes an useful integration time of 200 seconds. The 
pixel size A9 is 15".' With pixel array formats of 1024 x 
1024, observations cover a field of view of 4.27 degrees. 



use of the Fisher information matrix 



d 2 InL 

dpidpj 



(22) 



whose inverse provides the optimistic covariancc matrix for 
errors on the associated parameters (e.g., Tegmark et al. 
1997). In Eq. 22, L is the likelihood of observing data set 
x, in our case the angular power spectrum, given param- 
eters pi . . . p n involved in describing this data. Respect- 
ing the Cramer-Rao inequality (Kendall & Stuart 1969), 
no unbiased method can measure the ith parameter with 
standard deviation less than (F^) -1 / 2 if all other parame- 
ters are known exactly, and less than [(F -1 )^] 1 / 2 if other 
parameters are estimated from the data as well. 

For the present case involving C;'s of the IRB, we can 
write 



p y = ]T(Z + l/2)/ sky Tr 



idC t 1 dC l 



^P~ Cl 



d Pj 



(23) 



where / s k y = il/iir is the fraction of sky area covered by 
the imaging observations. Here, dCi/dpi denotes deriva- 
tives of the power spectra matrix, C^ A , with respect to 
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Fig. 7. — The expected errors on the reconstruction of star forma- 
tion rate based on IRB anisotropy angular power spectra measured 
with the proposed rocket-borne experiment. Here, we assume prior 
knowledge on the Pop III spectrum and the source bias; in the 
absence of knowledge on tgese parameters, one can constrain this 
combined quantity with the same fractional accuracy as shown in 
this figure. In general, the source spectrum may be considered to 
be known a priori due to the simplistic nature of first stars. The 
inner errors show the reconstruction assuming no foreground noise 
contribution, while the outer errors show a pessimistic case where 
we assume an additional shot-noise like power spectrum for the fore- 
ground noise, with a rms at the pixel scale of 10 nW m -2 sr — 1 . 



the parameter pi at each multipole and / s k y is the frac- 
tional sky coverage of the observations. This matrix can 
be written as 



'XX 



f £fnoise 



(24) 



where index i labels all contributions to the power spec- 
trum and (7 ; nolse is the instrumental noise contribution to 
the observations. Here, we assume that the noise contribu- 
tion between maps is uncorrelated. In terms of individual 
components of the power spectrum (ie. signal), we include 
contributions from the Pop III stars (which is considered 
to be the primary component of astrophysical interest) 
and the foreground galaxies. The clustering properties of 
zodiacal light is unknown and may be a major source of 
uncertainty in studies of early star formation with IRB 
fluctuations, but due to the unknown properties related 
to its spatial clustering, we do not consider its impact in 
this estimate. More over, the proposed multiwavelength 
data may allow us to reduce the confusion related to the 
zodiacal emission. 

In terms of parameters of interest, we study how well 
j\(z) can be studied with clustering information. We 
split this quantity in four redshift bins in the range of 10 
to 30 and extract the star-formation rate in those bins. 
Note that, in the strictest sense, what one extracts is 
(b(M))j\(z) which includes Pop III spectrum information 
and halo bias in addition to the star-formation rate. Here, 
for simplicity, we assume that both the spectrum and halo 
bias is known though one can constrain a combined quan- 
tity instead, if these latter two parameters are assumed 



to be unknown. Note the, in general, the spectrum, or at 
least the shape of the spectrum with an arbitrary normal- 
ization, can be assumed to be known apriori. 

We calculate the Fisher matrix to a maximum I of 5000 
such that the clustering of Pop III stars is in the linear 
regime (2-halo term) and is not shot-noise dominated (as 
shown in Fig. 3, we do not expect the Pop III shot-noise 
domination to be important until I ~ 10 4 ). In addition to 
the astrophysical parameters related to the star-formation 
rate, we also include cosmological parameters that define 
the linear matter power spectrum, mainly the power spec- 
trum normalization, tilt, Vl m h 2 and 0;,/i 2 . We set priors 
for these parameters as known from current WMAP re- 
sults and its extension related to information from other 
large scale structure surveys. This information, in gen- 
eral, leads to a power spectrum at accuracy at the level of 
few percent and is not a source of major concern for these 
studies. 

We summarize our results in Fig. 7, where we plot the 
star formation rate of the Pop III population and recon- 
structed errors on this rate based on clustering informa- 
tion for the parameters as tabulated in Table 1. We as- 
sume a total number of 512 2 pixels per band and consider 
two pointings with each pointing of 4.3 degrees field of 
view. At redshifts or order 10, where most of the Pop III 
contribution arises, the star formation rate is constructed 
with accuracies of order ~ 0.2 M sun yr" 1 Mpc~ 3 while 
at higher redshifts, ~ 30, the clustering information is no 
longer sensitive to the star formation rate due to the de- 
crease in fractional distribution of Pop III stars at such 
high redshifts. 

Note that we have assumed Pop III sources follow the 
optimistic model with a clustering angular power spectrum 
given by the upper curve shown in Fig. 4. If clustering were 
to be lower, then the reconstructed star formation errors 
increase. In such a scenario, observations with ASTRO-F 
can be used to investigate the existence of Pop III cluster- 
ing at arcminute scales to a lower flux level with a better 
removal of the foreground galaxies that may confuse the 
detection of Pop III sources. 

4. SUMMARY 

The recent results related to cosmic microwave back- 
ground (CMB) anisotropics suggest that the universe was 
reionized at a redshift around 17 with an optical depth for 
electron-scattering of 0.17 ± 0.04. Such an early reion- 
ization could arise through the ionizing radiation emitted 
by metal-free population III stars at redshifts of 15 and 
higher. We discuss the contribution to the infrared back- 
ground (IRB) surface brightness anisotropics from such a 
generation of early star formation. We have shown that 
the spatial clustering of these stars at tens of arcminute 
scales can potentially generate a significant contribution 
to the angular power spectrum of the IRB. 

Note that measurements of IRB spatial clustering al- 
ready indicate both an overall excess (e.g., Kashlinsky 
et al. 1996) and a specific excess at 100 arcmin scales 
(Matsumoto 2000, 2001) when compared to other angular 
scales. The models based on a contribution from Pop III 
stars indicate that one should expect a clustering excess 
at few tens of arcminutes related to the overall projection 
of the linear power spectrum at redshifts between 10 and 
30. While the amplitude of the clustering excess is model 
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dependent, the angular scale at which this clustering is ex- 
pected is more precise given that it is simply a reflection of 
the projected linear power spectrum at redshifts of 10 to 
30. The direct detection of the amplitude of clustering ex- 
cess, as a function of wavelength, can be used as a measure 
of the Pop III starformation rate, under the assumption of 
a spectrum and clustering bias, or the combined quan- 
tity involving the star-formation rate and bias weighted 
by the source flux. As an initial attempt in this direction, 
we have considered the astrophysical uses of a wide-field 
camera optimized for the detection of the first star signa- 
ture and have suggested that studies related to Pop III 
stars can be further improved with near-infrared obser- 
vations with the upcoming ASTRO-F mission. The best 
opportunity to study Pop III stars via IRB spatial fluctu- 
ations will come from a combined analysis of ASTRO-F 
and wide-field rocket-borne data in the same area on the 
sky. The higher resolution data from ASTRO-F can be 
used to clean low resolution, but significantly wider field 
of view, data from the rocket experiment, such that clus- 
tering studies can be extend even below our pessimistic 
level of the Pop III clustering considered. 
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